Femtometer Toroidal Structures in Nuclei 



J. L. ForestE and V. R. PandharipandeS 

Department of Physics, University of Illinois at Urbana- Champaign, Urbana, IL 61801 

Steven C. Piepeii and R. B. Wiringai 

Physics Division, Argonne National Laboratory, Argonne, IL 60439 

R. SchiavillaS 

CEBAF Theory Group, Newport News, VA 23606, 

and 

Department of Physics, Old Dominion University, Norfolk, VA 23529 

A. Arriag a0 

Centro de Fisica Nuclear da Universidade de Lisboa, Avenida Gama Pinto 2, 1699 Lisboa, 

Portugal 
(February 9, 2008) 



1 



Abstract 



The two- nucleon density distributions in states with isospin T = 0, spin 5=1 
and projection Mg=0 and ±1 are studied in 2 H, 3 ' 4 He, 6 ' 7 Li and 16 0. The 
equidensity surfaces for M$=0 distributions are found to be toroidal in shape, 
while those of Mg=±l have dumbbell shapes at large density. The dumb- 
bell shapes are generated by rotating tori. The toroidal shapes indicate that 
the tensor correlations have near maximal strength at r < 2 fm in all these 
nuclei. They provide new insights and simple explanations of the structure 
and electromagnetic form factors of the deuteron, the quasi-deuteron model, 
and the dp, dd and ad L=2 (D-w&ve) components in 3 He, 4 He and 6 Li. The 
toroidal distribution has a maximum-density diameter of ~1 fm and a half- 
maximum density thickness of ~0.9 fm. Many realistic models of nuclear 
forces predict these values, which are supported by the observed electromag- 
netic form factors of the deuteron, and also predicted by classical Skyrme 
effective Lagrangians, related to QCD in the limit of infinite colors. Due to 
the rather small size of this structure, it could have a revealing relation to 
certain aspects of QCD. Experiments to probe this structure and its effects 
in nuclei are suggested. Pair distribution functions in other T, S channels are 
also discussed; those in T, S = 1, 1 have anisotropics expected from one-pion 
exchange interactions. The tensor correlations in T, S = 0, 1 states are found 
to deplete the number of T, S = 1,0 pairs in nuclei and cause a reduction in 

nuclear binding energies via many-body effects. 
PACS numbers: 21.30.+y, 21.45,+v 
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I. INTRODUCTION 



Nuclear structure has been discussed mostly in the context of the liquid drop and shell 
models. These models have been extremely useful in explaining many observed nuclear 
properties. However, they are based on macroscopic concepts, and do not address the 
simplest nuclei, hydrogen and helium. Furthermore, recent (e, e'p) experiments [|I] have 
indicated that in heavier near-closed-shell nuclei, less than 70% of the nucleons are in the 
single-particle orbitals that would be fully occupied in the simple shell model. 

To obtain a more microscopic description of nuclear structure, we may regard the nucleus 
as a collection of interacting nucleons described by the Hamiltonian |2| 

^ = E-|^v l 2 + E^+ E (i.i) 

i i<j i<j<k 

The interactions u^- and V^k are not exactly known, but Vij is well constrained by the avail- 
able scattering data, and binding energies and theoretical considerations place important 
constraints on V^. The structure of the ground-state wave function \l/o a t small interparti- 
cle distances is influenced by the repulsive core and tensor parts of Vij. Most realistic models 
of contain these components and, for example, the Reid 0, Paris Q, Urbana [||] and 
new Argonne i>ig || models seem to predict similar structures. The three-nucleon interaction 
Vijk is much weaker than the v^. Due to a large cancellation between kinetic and two-body 
interaction energies, it has a significant effect on nuclear binding energies |2|] but its effect 
on the structure of \I/o is much less than that of the better known v^j. 

Due to the strong spin-isospin dependence of and Vijk it is difficult to solve the 
Schrodinger equation with the Hamiltonian (1.1). Only recently it has been possible to 
obtain accurate solutions for A < 7 nuclei J7]|| with the Green's Function Monte Carlo 
(GFMC) method. Accurate variational wave functions ty v , which contain less than 0.5% 
admixture of excited states are known for A = 3 and 4 nuclei ||. The available ty v for 
6 ' 7 Li [[Hj and 16 ||i~T| , |T2|1 are certainly not as accurate as those for A = 3 and 4, nevertheless 
they presumably contain most of the important structure of the \tV 
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In this paper we examine the short-range structure of 2 H, 3,4 He, 6,7 Li and 16 by cal- 
culating the two-nucleon density distributions in states with isospin T, spin S, and spin 
projection Ms- Since the deuteron has only two nucleons, its one- and two-body density 
distributions are trivially related. Variational wave functions *ff v and Monte Carlo methods 
are used for A > 2. 

The two-nucleon distributions in the T, S = 0,1 states have a strong dependence on 
the spin projection Ms- The equidensity surfaces, spanning the top three quarters of the 
density range in Ms=0 states, have toroidal shape. These tori are produced by the joint 
action of the repulsive core and tensor interactions. In contrast the equidensity shapes in 
the Ms=±l have dumbbell shapes, which have been studied earlier in the deuteron |]il||i"4"l . 
A brief description of the two-nucleon interaction in T, S = 0, 1 states is given in section 
II, and the density distribution of the deuteron is discussed in detail in section III, where 
we show that the dumbbell- shaped distributions in M$=±l states are produced by rotating 
tori. Commonly used models of Vij predict that the maximum density torus has a diameter 
of ~1 fm, and the half-maximum density torus has a thickness of ~0.9 fm. In section 
III we relate these dimensions of the toroidal distribution to the observed electromagnetic 
form factors of the deuteron. The structures are rather dense; current models predict the 
maximum one-body density of the torus inside the deuteron to be ~ 0.34 ± 0.02 fm -3 , i.e., 
approximately twice nuclear matter density. 

The two-nucleon T,S = 0,1, Ms = 0, ±1 density distributions in 3 ' 4 He, 6,7 Li and 16 
are compared with those of the deuteron in section IV. The distributions for r < 2 fm differ 
only by a single scale factor. They indicate that in the T = state, the tensor correlations 
have near maximal strength for r < 2 fm in all these nuclei. The scale factor is identified as 
the Levinger-Bethe quasi-deuteron number, and its value is compared to the ratios of total 
photon (i? 7 = 80 to 120 MeV) and pion (E 7T + ~115 MeV) absorption cross sections. 

In order to study the nature of many-body structures induced by these compact two- 
body structures we study the dp, dd and ad overlaps with the ty v of 3 He, 4 He and 6 Li in 
section V. These depend strongly on the spin projection Md of the deuteron and indicate the 
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presence of anisotropic structures in all these nuclei. Experiments to probe these structures 
are suggested. 

Pair distribution functions in other T, S, Ms states are discussed in section VI. Those in 
T,S — 1,1, Ms = 0, ±1 states are anisotropic as expected from the pion-exchange tensor 
force. We also find that the number of T, S = 1,0 pairs in a nucleus is reduced due to 
many-body effects involving the strong T, S = 0, 1 tensor correlations. This reduction gives 
a significant contribution to the saturation of nuclear binding energies. 

The Skyrme field theory [|T5 1 , related to QCD in the limit of large number of colors 
N c — > oo, has predicted toroidal shapes for the deuteron |16H18| in the classical limit. 



Density distributions of the ground states with 3 to 6 baryons have also been calculated [19 



in this limit. In the last section, VII, we summarize our results, obtained with conventional 
nuclear many-body theory, and indicate their relation to those of the Skyrme field theory in 
the classical limit, and of the constituent quark model. 



II. THE TWO-NUCLEON INTERACTION IN THE T, S = 0, 1 STATE 

Nuclear forces are not yet quantitatively understood from QCD. However, many realistic 
models have been constructed by fitting the available two-nucleon scattering data. The 
shape of the short-range structures in the T, S = 0, 1 state appears to be relatively model- 
independent. The interaction t>o,i m the T, S = 0, 1 state in Reid, Urbana and Argonne 
models can be expressed as 

vo,i = «g fl (r) + «S fl (r)5 y + < x (r)L • S + < x (r)L 2 + < 2 (r)(L • S) 2 , (2.1) 

while a V 2 operator is used instead of L 2 in the Paris potential. The structures are formed 
mostly by the static part of the interaction: 

<f =<i(r)+<!(r)% (2.2) 

It is instructive to study the expectation values of v^f in eigenstates of the position operator 
r with spin projections Ms=0 and Mg=±l. These depend upon r and 8, the polar angle of 
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r with respect to the spin-quantization axis z, and are given by 

(M s =0\v s f(r)\M s =0) = v c 0tl {r) - 4< 1 (r)P 2 (cos 9), (2.3) 
(M s = ± 1\v$*(t)\M s = ± 1) = v c 0A (r) + 2v* il (r)P 2 (cos0). (2.4) 

The Ms=0 expectation value of fjff* has the largest variation with 9 as illustrated in Fig. [TJ. 
The static potential has a repulsive core; outside the core it is very attractive for 9=tt/2 
and repulsive for #=0 and tt. Therefore, in this state the np pairs form a toroidal density 
distribution in the xy plane. The potential in the M$=±l states is attractive for #=0 and 
7r, and equal to that for M$=0, 9=n/2, while it is repulsive for 9=ir/2 and half way between 
the Ms=0, 9=0 and 7r/2 potentials. Thus the Mg=±l potential has two distinct minima 
separated by a barrier, and therefore the density distributions have a dumbbell shape in this 
state. 

At r >1.5 fm, the v^f is dominated by the one-pion-exchange potential, while at smaller 
r it has a significant model dependence. Much of this model dependence is cancelled by the 
differences in the momentum-dependent terms in the models. In particular, the deuteron 
wave functions calculated from these potentials have much smaller model dependence. These 
are commonly written as 

*f d (r) = R (r) y^(f) + R 2 (r) 3#?(r), (2.5) 

where Ro{=u/r) and R<2(=w/r) are the S- and D-state radial wave functions and yffsj are 
the spin-angle functions. The Ro and R2 calculated from the different potential models are 
shown in Fig. The short-range structures are related to the Rq and R2 functions, and 
therefore we expect them to be fairly model- independent. 

However, the "full" Bonn potential |2(J offers an exception. The Rq and R2 predicted 



by this one-boson-exchange model of the A^A^ interaction are similar to other predictions at 
larger r, but they have an additional sharp structure with a range of ~0.2 fm. We will not 
consider the possibility of such an additional structure in this work. 
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III. THE DEUTERON 



The short-range structure of the deuteron is most obvious in its density distribution 
p^f d (r', 9) which depends upon the projection M d of the total deuteron angular momentum, 
the distance r' from the deuteron center of mass, and the polar angle 9 of r'; it is independent 
of the azimuthal angle (p. Note that the interparticle distance r = 2r', and the standard 



normalizations, 



2 



r dr 



o 



R 2 o(r) + RUr)] =1, (3.1) 
|dVpfV) = 2, (3.2) 



are used in this work. 



The p d Id (r') is given by I6^ d d (2r')#f d (2r'), where the factor 16 comes from the differ- 
ence in normalizations (3.1) and (3.2). A simple calculation using Eq. (2.5) yields 

Pd (r') = -[C (2r')-2C 2 (2r')P 2 (cos9)\ , (3.3) 

7T 

PdV) = - [C o (2r0 + C 2 (2r')P 2 (cos£)] , (3.4) 

7T 

with 

C (r) = R 2 (r) + Rl(r) , (3.5) 
C 2 (r) = V2Ro(r)R 2 (r)-^Rl(r) . (3.6) 

The interesting structure of these density distributions is shown in Figs. |3|-|| Fig. ^| shows 
p'd d i r ') along 9 = and 9 = tt/2 directions, noting that 

p d (r',9=n/2)= P r(r',9=0). (3.7) 

The above densities are the largest while p°(r', #=0) is the smallest as expected from the 
properties of v^f discussed in Sec. II. The small value of the ratio p d (r',9 = 0)/p d (r',9 = 
7r/2) indicates that the deuteron has near maximal tensor correlation at distance r' < 1 fm 
or equivalently at r < 2 fm. This ratio is ~0 for maximal tensor correlations. 



Figures £| and || show the density distributions predicted by the Argonne Vi$ model in the 
x'z' plane. The maximum value of p d is fairly model-independent (Fig. |3]) and large (~0.35 
fm -3 ). The maxima of p° d (Fig. |j) form a ring with a diameter of ~1 fm, denoted by d, in 
the x'y' plane, while the p^ 1 has two equal maxima on the z'-axis (Fig. ^|) at z'=±d/2. 

The three-dimensional distributions p^ d (r') can be obtained by rotating the distributions 
shown in Figs. f| and || about the z'-axis. They are represented by equidensity surfaces shown 
in Fig. |] for p d =0.24 and 0.08 fm -3 ; all four sections are drawn to the same scale. The 
maximum value of p d (8=0) is ~0.05 fm -3 (Fig. |||). Therefore the equidensity surfaces for p d 
having p d >0.05 fm -3 cannot intersect the z'-axis, and thus have toroidal shapes shown in 
Figs. |6|B and |6p. The central hole in these tori is due to the repulsive core in Vq^, and their 
angular confinement is due to the tensor force. In absence of the tensor force, -R2O") = 0, 
the p d = p d l , and the equidensity surfaces are concentric spheres. 

The maximum value of p d zl {9= r K /2) is ~0.19 fm -3 as can be seen from Fig. |3|. Therefore 
the equidensity surfaces of p^ 1 for p d >0.19 fm -3 can not cross the 9=tt/2 plane; they have 
two disconnected parts forming a dumbbell as shown in Fig. ^]A. At smaller values of p d we 
also obtain two equidensity p^ 1 surfaces (Fig. |5]C), consisting of an inner surface due to the 
repulsive core enclosed by an outer. At very small p d (<0.05 fm -3 ) the equidensity surfaces 
of p° d also have disconnected inner and outer parts, neither close to spherical in shape. 

The toroidal shape of the M d =0 equidensity surfaces is more compact and persists down 
to smaller p d , or equivalently to larger values of r', as can be seen from Figs. |||| In the 
classical Skyrmion field theory only this shape is obtained for the distribution of baryon 
density in the ground state for two baryons [|1^]. The deuteron can be considered to be 
more deformed in the M d =0 state. For example, the expectation values of the quadrupole 
operator Q(r') = 2 r /2 P 2 (cos#) obey the relation 

/ p°(r') Q(r') dY = -2 J pf(r') Q(r') dY. (3.8) 

It is rather simple to obtain the p^ :1 (r / ) from the p d (r'). We rotate the p°(r') about the 
y'-axis by an angle of tt/2 so that the toroidal ring is in y'z' plane with x' as the symmetry 
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axis. This places the deuteron in the superposition of M^=±l states. The Mrf=±l states 
are obtained by spinning the rotated toroid about the z'-axis, and the p^ l (r') is just the 
average value of the density of the spinning toroid, i.e., 

pf(r\d) = — r p° d (r', cos" 1 (sin0 cos <f>)) d0. (3.9) 
2% Jo 

The L=l, M=±l states of the harmonic oscillator, given by (f){r)=e~ vr (x±iy), are obtained 
in the same way from the (f)(r)=e~ ur2 z, L=l, M=0 state. Therefore it is tempting to consider 
the toroidal shape of p° d as the basic shape of the deuteron. The expectation value of the 
current operator is zero in the M d =0 state, therefore one may regard that as the "static" 
state of the deuteron. Note that the toroidal shapes cannot be obtained by rotating the 
dumbbell by tt/2 about the y'-axis and spinning it about the z'-axis. This gives p d (8=0, r') = 
p d l (6=ir/2,r') which is not true. The dumbbell- or cigar-shaped density distribution of the 
deuteron in the M d =±l state has been studied earlier fl3|]l4|| . Unfortunately the toroidal 
distribution of the M d =0 state was not studied, and its similarity with the predictions of 
Skyrmion field theory was not noticed. 

The deuteron electromagnetic structure functions A(q) and B(q), and the tensor polar- 
ization T 2 o(q) in elastic electron-deuteron scattering have been extensively studied exper- 
imentally [^l| -|30[| and theoretically j3lH34f. They are usually calculated from the S- and 



D-wave functions Rq and R 2 obtained from realistic interactions, by including in the nuclear 
electromagnetic current, in addition to the dominant impulse approximation (IA) operators, 
relativistic corrections and two-body meson-exchange contributions [pT] , p3[ . More recently, 



calculations of these observables based on quasipotential reductions of the Bethe-Salpeter 
equation and one-boson-exchange interaction models, constrained to fit nucleon-nucleon 
data, have also been carried out p^j3^] . The theoretical predictions for the structure func- 



tions based on both the nonrelativistic and relativistic approaches are in good agreement 
with data. Our interest here is not to improve upon the present theoretical predictions, 
but to relate the values of the minima and maxima of T 2 o{q) and B(q) to the size of the 
toroidal structure in the deuteron. The q- values of these extrema may be shifted by ~10% 
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by corrections to the impulse approximation used in the following pedagogical analysis. 
The charge form factors, defined as 

1 r n M dfVW^'j3/ 



Fc,mM) = 2/ ^ V)e^dV, (3.10) 

are shown in Fig. [7[ At large q the Fc } \{q) gets most of its contribution from the two 
peaks of p^(r') (Fig. [|) at z'=±d/2. The Fourier transform of the sum of two (^-functions 
at z'=±d/2 is given by cos(qd/2) with zeros at qd=ir, 3tt, ■ ■ ■. These zeros are due to the 
cancellation of the contribution from the two peaks, and they persist even when the peaks 
have a finite width. The first two minima of F^^q), obtained from the Argonne v\$ p^(r'), 
occur at gi=3.6 and ^2=12. 6 fm -1 . The effective values di estimated from the minima g^, 
using di=(2i — l)7r/gj, are 0.87 and 0.75 fm for 2. These values are smaller than the 
diameter d=l fm because the dumbbell ends (Fig. ^A) are not spherical. Nevertheless the 
minima of x seem to be primarily determined by the diameter d of the maximum-density 
torus. 

The Fourier transforms of a disc of thickness t, with q perpendicular to the disc, are 
proportional to sin(qt/2)/(qt/2) irrespective of the shape of the disc. These have zeros at 
qt=2ir, 4ir, ■ ■ •, which may be used to obtain the thickness t. The first two zeros of Fc t o(q) 
at g=9.2 and 19.5 fm" 1 (Fig. 0) give values 0.68 and 0.64 fm for the effective thickness of 
the torus. The maximum thickness along the z-axis of the calculated equidensity surface at 
half maximum density is 0.88 fm. 

The T 2 o{q) form factor of the deuteron has small magnetic contributions which depend 
upon the electron scattering angle 9. The extrema of T 2 o are not significantly affected by 
this magnetic contribution as can be seen from Fig. [8], and ignoring it we obtain a rather 
simple equation: 

rF*Jq) -F^Jq) 

Its minima occur when F ( S 1 (g)=0, while the maxima have F^, (q)=0. These minima and 
maxima correspond to those values of q at which the recoiling deuteron has only Md=0 
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or Mrf=±l, respectively. The first minimum of T 2 q is experimentally known to occur at 
q ~3.5±0.5 fm -1 in agreement with the value d ~1 fm predicted by realistic potentials. The 
first maximum of T 2 $ has not yet been experimentally located; it provides a measure of the 
thickness t. 

In magnetic elastic scattering the deuteron spin projection in the q direction changes 
by ±1 since the photon has M 7 =±l. Thus the magnetic form factor Fm{q) is a transition 
form factor. It has convection current and spin-flip terms [j3T| , ^3|| of which the latter is 
dominant. The Fm{q) calculated with and without the convection current term are not too 
different (Fig. |]). The spin flip part of Fm{q) can be obtained from the transition density 
Ar(r'): 



F s J(q) = (n P + Hn) I e^'p^rOdV 



Ptr (r') = - [R 2 (2r>) - -Rl(2r') - - 



V2R (2r')R 2 (2r') + R 2 2 (2r')} P 2 (cos0) 



(3.12) 
(3.13) 



This transition density is shown in Fig. |H| it is dominated by the toroidal p°(r'), and its 
effective thickness along z' axis, obtained from the zeros of Fm{q) (Fig. 0), is ~0.85 fm. The 
minimum of Fm{q) is observed |26| at q ~ 7 fm -1 , supporting the theoretical estimates of t. 
The deuteron wave function in momentum space is defined as 



(2tt) 3 / 2 

Ro(k) y^(k) + R 2 (k) y™f(k) 



(3.14) 



with 



R L (k) = i L J2/7T / dr r 2 j L (kr)R L (r) . 



(3.15) 



The momentum distributions p^ d (k), given by ^/^ d (k)^f d (k), are then easily obtained as 

1 



(k) 



47T 
1 



C (A;)-2C , 2 (A;)P 2 (cos^) 



(k) = — C (k) + C 2 (k)P 2 (cos8 k ) 

47T L 



(3.16) 
(3.17) 



where 9 k is the angle between k and z-axis, and the Ci(k) are defined as in Eqs. (3.5)-(3.6) 
with Rl{t) replaced by Ri{k). Note that the pd are normalized such that 
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/d 3 kpf d (k) = l • (3.18) 
The momentum distributions p®(k,9k) and p^ x {k^9k) for Oy. = and tt/2 are shown in 



Fig. [LJ. Note that 

~p° d (k,e k =n/2) = ~pf(k,e k =0) . (3.19) 

The zeros of f>d(k, 9 k =0) and p%k, 0^=0) occur at k ~ 1.5 fm -1 and 5.2 fm _1 , respectively, 
and are related to the spatial dimensions of the torus. In naive estimates these minima 
occur at 7c/2d and n/t, respectively. Thus measuring the positions of the zeros in these 
momentum distributions would provide an independent estimate of the spatial dimensions 
of the toroidal structure in the deuteron. This information would be complementary to that 
yielded by elastic form factors measurements. 

The momentum distributions p^ Id (k,9k) could in principle be measured by (e,e'p) scat- 
tering on polarized deuterons. In the one-photon exchange approximation the d(e, e'p)n 

cross section, in the laboratory frame, is generally expressed as 
d 5 a M d 



dE' e dn' e dn p 



T7 _ E P E n 



a MoU PpEpR Te l {v L R L + v t Rt + v lt Rlt + v T tRtt) , (3.20) 

EpP r , 



m 2 



Pp • Pn 



(3.21) 



E n p p 

where is the target spin projection, E' e and Q' e are the energy and solid angle of the final 
electron, and fl p is the solid angle of the ejected proton. The coefficients v a are defined in 
terms of the electron variables, while the structure functions R a involve matrix elements 

(n + p; p„p p , M n M p | L)T (q) | d, M d ) , (3.22) 

of the charge (Ol) and current (Ot) operators between the initial deuteron and final n + p 
states. The neutron momentum is p n =q — p p , q is the momentum transfer, p m =-p n is 
the missing momentum, and M p and M n are the proton and neutron spin projections, 
respectively |3"3] . The cross section for unpolarized deuterons, 

d 5 a 1 ^ d 5 a Md 



dE' e dQ' e dn p 3 Md t^±i dE' e dQ' e dn p 
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has been measured up to p m ~500 MeV/c, and there is good agreement between theory and 



experiment |35 



In plane- wave-impulse-approximation (PWIA), obtained by neglecting interaction effects 
in the final n + p states as well as relativistic corrections and two-body terms in the charge 
and current operators, the M rf -dependent cross section is proportional to 

d 5 a M " 



cxpf*( Pm ) . (3.24) 



dE' e dVt' e dVt p 



Using, for example, tensor polarized deuterium, it should be possible to measure the differ- 
ence between p°(p m ) and (p m ), and therefore empirically determine the positions of the 
minima in these momentum distributions. Clearly, such an analysis is justified if the PWIA 
is valid. This assumption has been tested by carrying out the full and PWIA calculations of 
the d 5 a Md / dE' e dQ' e dQ p in parallel kinematics with q fixed at 500 MeV/c, u in the range 290- 
390 MeV, and the electron scattering angle 8 e = 10°. The results, shown in Fig. [12|, indicate 
that, while FSI, two-body current, and relativistic corrections are not entirely negligible, at 
least in the kinematical region which has been studied here, their effect is small compared to 
the difference between the cross sections for Md—0 and ±1. We therefore conclude that the 
results of such an experiment could be used to empirically study the diameter and thickness 
of the torus. One might argue that this information could be more easily obtained from 
elastic form factors measurements, as discussed above. However, it should be realized that, 
in contrast to the d(e, e')d data, the double-coincidence data would allow us to ascertain to 
what extent this toroidal structure is due to nucleonic degrees of freedom. 

IV. THE TWO-NUCLEON DENSITY DISTRIBUTIONS IN NUCLEI 

The two-nucleon density-distributions in T, S, M T , M s two-nucleon states are defined 

as 

P^ s ' Ms (v) = — ±— f (^| E^T.^Mt.Ms)^) . (4.1) 

ZJ + 1 Mj=-J i<j 
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Here |\l/ Mj ) denotes the ground state of the nucleus with total angular momentum J and 
projection Mj, and Py (r, T, S, M T , M s ) projects out the specific two-nucleon state with 
Ti — Tj = r. For N = Z nuclei, the wave functions used in this study are symmetric under 
exchange of neutrons and protons; hence Pt t s Ms { t ) is independent of Mt- For 3 He, we have 

p5 Ms = 0; P tT s = ^ s , (4-2) 

while for larger N ^ Z nuclei, the Mt dependence is nontrivial. In the following we discuss 
Pt,si the average over Mt of PtJs' Ms ■ The p^s is normalized such that 

E (2T+1) / p A T %r) d 3 r = \a(A-1) , (4.3) 

T,S,Ms 

which is the number of pairs in the nucleus. It is a function of r and 9 independent of the 
azimuthal angle 0. 

It can be verified from Eqs. (3.3-3.6) and (4.1) that in the deuteron 

PoS( r ) = ^x^V = r/2) . (4 . 4) 

Note that the spin-dependent two-body density on the left (Ms = M) is an average over 
projections M^ in the deuteron, while the polarized one-body density on the right (M^ = M) 
has been summed over spins. The p^ s in 3,4 He, 6 ' 7 Li, and 16 O have been calculated from 
variational wave functions using Monte Carlo techniques. For the A < 7 nuclei, these 
wave functions minimize the expectation value of a Hamiltonian consisting of the Argonne 
vig two-nucleon and Urbana IX three-nucleon potentials (for A = 6,7 the minimization 
is constrained by the experimental rms radii); a detailed description of the form of the 
wave functions is given in Refs. fl]10||. The wave function for 16 O was obtained from the 
variationally best wave function by slightly increasing the radius of the single-particle part 
of the wave function so as to reproduce the experimental rms radius of 16 O. The details of 
the 16 O wave function will be published elsewhere ||12|| . A cluster-expansion including up 



to four-body clusters with Monte Carlo integration [II] was used to compute the two-body 
densities in 16 O. 
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To reduce statistical fluctuations in the calculation, we write 



P*Stt= E ^§»P L (cos0) , (4.5) 



L=0,2 



and directly compute the A(r) as 

= 27V1 ^ x 

2 / dR* M '(R)t £ ^5(r - r y ) P L (r y • 5) P^T, 5, M s ) *^(R) , (4.6) 

where R represents the coordinates ri, . . . , r^. Because of the average over the total spin of 
the nucleus, the A^ s s L are zero for L > 2, and 

A M S =0 _ A M S =±1 i. ? x 

A M S =0 __n A M s =±l (A o\ 

For the deuteron, the A^ L are related to the C L of Eqs. (3.3-3.6) by A$f fl = C /48; A\ A2 = 
C 2 /48. 

The shapes of p Q f (r, 9) are very similar at r < 2 fm in all the nuclei considered. In order 
to study the evolution of p^f with A we divided the p i(r, @)a by the ratio R A d defined as 



Max (ptfae: 



The densities so normalized are compared in Fig. 13, and the values of Raci are listed in 



Table I. Fig. H shows p° 0A (r, 9 = 0) A /R Ad , pg(r, 9 = n/2) A /R Ad , and p° x (r, = n/2) A /R Ad 
for 2 H, 4 He, and 16 0. Note that p[] 1 (r, 6 1 = ir/2) = Po\(r,9 = 0) by virtue of equations 
(4.4), (4.6) and (4.7) in all nuclei. After normalization by R A d, the various densities for 
3 He lie between those of 2 H and 4 He, while those for 6,7 Li are in between the 4 He and 16 
results. It is obvious from Fig. O that the equidensity surfaces of the two-body density p^ff 



are very similar to those of the deuteron density shown in figures at r < 2 fm (r' < 1 
fm). At r < 2 fm the ratio p[] 1 (r, 6 1 = 0)/p[] 1 (r, 9 = tt/2) is very small, indicating that the 
tensor correlations have near maximal strength in all the nuclei considered. In 16 O the p f 
becomes approximately independent of Ms only for r > 3 fm. 
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Bethe and Levinger suggested in 1950 [36| that at small distances the relative T, S 



0,1 neutron-proton wave function in a nucleus is likely to be similar to that in the deuteron. 
We find that this is a good approximation. The expectation value of any short-ranged two- 
body operator that is large only in the T, S — 0, 1 state scales as R,Ad- In Table I we list 
values of the ratios of the calculated expectation values of the one-pion exchange part of 
the Argonne vi$ potential, the observed low-energy (118 MeV for 3 He |37] and 4 He [[38[| , and 



115 MeV for 16 O |39|]) pion absorption cross sections and the average value of the observed 
photon absorption cross sections in the range E 1 = 80 to 120 MeV. All these processes are 
dominated by the T, S — 0, 1 pairs, and seem to scale as RAd- 

While comparing these ratios in detail it should be realized that (v n ) in nuclei has a 
relatively small contribution from T, S ^ 0, 1 states, absent in the deuteron, which makes 
(v n ) a/ (v-n) d slightly larger than RAd- The (v n )d = -21.3 MeV for the Argonne t>i 8 model, 
and it accounts for most of the deuteron potential energy, (v)d = -22.1 MeV. Also in larger 
nuclei, the {v^)a gives a large fraction of the total two-body interaction energy ||11|| . Direct 
comparison of the ratio of pion absorption cross sections with RAd may not be strictly valid. 
The scattering and absorption of pions by spectator nucleons, absent in the deuteron, is 
expected to reduce the ratio cr^ b A /a1 b d1 while three-body and higher absorption mechanisms, 
also absent in the deuteron, will increase it. After correcting for final state interactions of 
the two outgoing protons, the two-body ,pp) part is estimated to account for ~ 76% of 
the total absorption cross section for 115 MeV tt + by 16 3~9|| . In 3 He ~ 20% of the 118 



MeV 7r + absorption cross section has three-body character [37]], however a part of this 20% 
must be due to initial and final state interactions. 

Results of Mainz experiments |4(J on 7 Li and 16 O are used to calculate the average value 



of cr 7 in the energy interval E y = 80 to 120 MeV. The cr^, d averaged over the same energy 
interval is ~ 0.072 mb The only available data for 3 He in this energy range are from 
the experiments done in the 1960's and 1970's [|3|. The average cross section of the 



two-body photodisintegration of He, in the energy range 80-120 MeV, is ~0.03 mb |H3 



and that for the three-body process is ~0.10 mb [)Eql , giving total cross section of ~0.13 mb. 
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The average cross section for total absorption of photons by 4 He in the same energy range 
is crudely estimated from Fig. [I] in ref. to be ~0.3 mb. 

The total number of pairs with given T, S in nuclei can be computed as 

N T,S = £ ( 2T + !) 21T I r2dr dc ° s6 PT,s( r > d ) A 
M S 

= (2T + l)(2S + 1) 4tt / r 2 rfr A° j50 (r) ; (4.10) 

the values for T, S — 0, 1 and the corresponding naive independent-particle model values 
are also shown in Table I. We see that the correlations induced by the potentials do not 
significantly change the from their independent-particle (IP) values; however, as will be 
discussed later, this is not true for T=l pairs. For few-body nuclei, Ra<i is significantly larger 
than Nq X) however, in 16 O N$ x has a large contribution from pairs with large r and Racl 
is smaller than Nq X . The calculated value of R A d for 16 is much smaller than Levinger's 
estimate R Ad ~ 8NZ/A |6|. 

The pd I=o (r,0 = n/2) has its half-maximum value at r ~1.8 fm (Fig. |13|). If we identify 
the region with r < 1.8 fm as the "quasi-deuteron" , then the probability that the np pair in 
a deuteron is in the quasi-deuteron region is ~0.25, and the number of quasi-deuterons in a 
nucleus is ~ In the past, however, R^ d itself has been interpreted as the number of 

quasi-deuterons in the nucleus. 



V. TWO-CLUSTER DISTRIBUTION FUNCTIONS 

The strong spin-dependent anisotropics of the two-nucleon densities suggest that three- 
nucleon and higher distribution functions in nuclei could also be anisotropic. A general 
study of these higher distributions is beyond the scope of this work; however the two-cluster 
distribution functions dp in 3 He, dd in 4 He, and ad in 6 Li are simple to study with the Monte 
Carlo method |^5| . They provide some information on the higher distribution functions, and 



may be relatively accessible by (e, e'd) and (e, e'p) experiments. 
The two-cluster overlap function can be written generally as 
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A ab (M a , M b , Mj, r ab ) = (A^ A J a K b , r ab \* Mj ) (5.1) 

= J2 (LM L SMs\JMj)(J a M a J b M b \SM s )R L (r ab )Y LML {r ab ) , 

LM L SM S 

where r ab is the is the relative coordinate between the centers of mass of the two clusters and 
A is an antisymmetrization operator for the two-cluster state. The Rii^ab) radial functions 
can be evaluated from 

R L (r ab )= E (J a M a J b M b \SM s )(LM L SM s \JMj) 

M a M b M L M s 

f rfR[^^(R a )^^(R b )]ty; ML (f afc ) (5(r ~ ^V ^R) , (5.2) 

J r ab 

where R a (&) represents the coordinates of particles in cluster a{b). We note that in PWIA the 
(e,e'a)b cross section is proportional to the momentum distribution \A ab (M a , M b , Mj, k)| 2 
obtained from the Fourier transform of the overlap function A ab (M a , M b , Mj, r ab ). 

In the present work, the integrations have been made with Monte Carlo techniques akin 
to Ref. ||45|| , but with some improvements. Configurations are sampled with the weight 



function l 1 !/*^ 7 ! 2 containing the full variational wave function. In Ref. only a single term 
in the antisymmetric product in Eq. (5.2) is calculated. The efficiency of the Monte Carlo 
sampling has been improved by evaluating all possible partitions of the nucleus into clusters 
a and b at each configuration R. We also use a much larger sample of configurations than 
in the previous calculations. 

We can define a two-cluster wave function, in analogy with the deuteron wave function 
of Eq. (2.5), using the radial overlap functions 

*S , (r,*) = EW3ffi(^) 

LS 

= E A ab (M a ,M b ,Mj,r ab ) X Ma X Mb , (5.3) 

M a M b 

where \ Ma and \ Mb denote spin states J a M a and J b M b of a and b. For the cases ab = dp, 
dd, and ad there are both S- and D-wave states in the two-cluster wave function. In these 
cases the well known D 2 parameter can be defined by means of the Ro(r) and R2{r) radial 
functions [46[ : 



,o6 



/ R2{r ab )r 4 ab dr ab 



(5.4) 



15 J R (r ab )rl b dr ab ' 



Although in the present paper we emphasize the short-range structure of nuclei, it is also 
interesting to study the asymptotic behavior of the overlap radial functions. Of particular 
interest is the asymptotic D/S ratio r\ ab = C^/C^, where C and C 2 are the asymptotic 
normalization constants of Ro(r) and Ri{r) respectively: 



Here hi is the spherical Hankel function of first kind and a ab is the wave number associated 
with the separation energy of the nucleus into clusters a and b. We must point out that the 
present variational method, as well as the GFMC method, determine the wave functions by 
energy minimization, to which long-range configurations contribute very little. Therefore 
these methods are not very sensitive to the asymptotic part of the wave functions, and 
consequently our values for r\ ab should be considered only as estimates. 

The two-cluster density distribution for a given set of spin projections is defined as: 



In each of the cases studied here, it exhibits spin-dependent spatial anisotropies which are 
easily understood in terms of the toroidal or dumbbell structure of the polarized deuteron. 
The density is enhanced in the direction corresponding to the most efficient or compact 
placement of the deuteron with the remaining cluster, and reduced in those directions that 
would lead to very extended structures. 

Finally, we are also interested in the total normalizations N L of the S- and D-wave 
two-cluster distributions: 



These quantities can be related to spectroscopic factors and give the total S- and D-state 
fractions. All the results presented here are obtained from the Argonne v\$ two-nucleon and 
Urbana model IX three-nucleon interactions. 



RL{r a b) = lim -i L Cl b h L (ia ab r ab ). 



(5.5) 




(5.6) 




(5.7) 
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A. dp Distribution in 3 He 



The calculated Ro(r dp ) and ^(r^) are shown in Fig. |TJ|; the R2{r dp ) is negative and 
smaller in magnitude than the R 2 in deuteron (Fig. ^). The D 2 P value obtained with these 
radial functions is —0.15 ± 0.01 fm 2 , a little smaller than experimental estimates, ranging 
from -0.20 ± 0.04 to -0.25 ± 0.04, obtained through DWBA analysis of (ci, 3 He) transfer 



reactions . In Fig. 14 we also show our asymptotic fit to the S- and D-waves. The result is 
i] dp = —0.035, somewhat smaller than the Faddeev result |47j], —0.043 ±0.001. Experimental 
estimates, also obtained through DWBA analysis of (<i, 3 He) transfer reactions, range from 



-0.042 ±0.007 to -0.035 ±0.006 46 



The total normalizations are Nq P = 1.31 and N$ p = 0.022. The sum, 1.33, can be 
interpreted as the number of deuterons in 3 He | 4"5| . It is less than 1.49 (Table I), the number 
of T, S = 0, 1 pairs, because the pairs are not always in the deuteron state. It is also smaller 



than RAd = 2.0 inferred from short range distribution functions (Fig. 13 and Table I). This 
is probably because 3 He is more compact than the deuteron. 

The p dp ' 2 and p d 2 ' 2 are shown in Figs. |T5] and When M d = 0, the M p = ±~ 
proton is preferentially along the z-axis; in contrast, when M d = 1, the M p = — | proton 
is more likely to be in the xy plane. In the first density distribution the S- and D- wave 
amplitudes interfere constructively, to enhance the probability of finding the proton along 
the z-axis, whereas in the latter the interference is constructive on the transverse xy-pl&ne. 
Consequently _R and R 2 have opposite signs and and r] dp are both negative. The spin- 
dependent dp anisotropies are favored by both tensor and central forces, and lead to more 
compact three-body states. 

The momentum distribution of dp clusters in Mj = | 3 He is shown in Fig. [17| for M d , 
M p = 0, I and 1, — | for momenta parallel and transverse to the z-axis. In PWIA the 
3 He(e, e'p)d cross-section is directly related to these momentum distributions. A large spin 
dependence of the missing-momentum distribution for protons ejected in parallel kinematics 
is predicted for q parallel to z. The minimum for the momentum distribution along z-axis 
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for M p = +| occurs at ~1.4 fm _1 , while that for M p = — | is at ~2.4 fm _1 . Thus the spin 
asymmetry, (n-f — Ji^)/ (nj + of the protons ejected from polarized 3 He changes from ~ 
-1 to +1 as the missing momentum varies from ~1.4 to 2.4 fm -1 in parallel kinematics and 
PWIA. The dp momentum distribution in unpolarized 3 He has been studied at Saclay |48|] up 
to ~2.5 fm _1 . The observed distribution is generally smaller than the PWIA prediction |35 



indicating attenuation due to FSI. However, a part of the FSI attenuation will cancel in the 



asymmetry, and moreover, it is now possible to perform continuum Faddeev calculations f49 
including FSI. 



B. dd Distribution in 4 He 

The calculated Ro(r dd ) and R2{rdd) are shown in Fig. [Tj| The D dd value obtained with 



these radial functions is — 0.12±0.01 fm 2 . In Fig. [TJ]we also show our asymptotic fit to the S- 
and D-waves. The result is rjdd = —0.091. The integrals of these functions yield N^ d = 0.98 
and N dd = 0.024. The number of deuterons present is greater than twice the sum of these 
quantities, 2.0, when one allows for the additional presence of d + p + n configurations. The 
p dd '° and p^ '° are large in 4 He and their anisotropics, induced by the tensor interaction 
and the shapes of deuterons, are shown in Fig. [H| The p d '^° is largest when r dd is along 
the z-axis, i.e., when the deuterons are in the toroidal shape and have a common axis. It is 
smallest when r dd is transverse (two tori side by side) and equal to that for p dd 1,0 with r dd 
parallel to z (two dumbbells in a line). The latter distribution is of intermediate strength 
when r dd is transverse (two dumbbells side by side). Again in the first (second) density 
distribution the S- and D-wave amplitudes interfere constructively (destructively) along the 
z-axis. Therefore Ro and R2 have opposite signs and D dd and r\ dd are both negative. 

The momentum distributions are also anisotropic (Fig. |2(|). In particular the p dd '°(kz) 
has a dip at k ~1.7 fm -1 that is absent in the p dd 1,0 (kz). It may be possible to study 
these with (e, e'd) reactions. The unpolarized 4 He(e, e'd)d reaction has been studied at 
NIKHEF |50|| . The observed cross sections are much smaller than estimates using the dd 
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momentum distribution and either PWIA or DWIA. 



C. ad Distribution in 6 Li 



The calculated Ro(r ad ) and R2(r a d) are shown in Fig. BTl The Ro(r ad ) and R^fad) both 



exhibit nodes at short distances and have opposite signs almost everywhere. This nodal 
structure has been predicted ina + d and a + p + n cluster models, but not always with 
the correct relative sign The asymptotic behavior is correlated with the quadrupole 

moment Q; obtaining the experimental value of —0.08 fm 2 has been a long-standing problem 
in a + p + n cluster models. The variational wave function used here gives Q = —0.8 ± 0.2 
fm 2 , i.e., the correct sign but far too large in magnitude. Small changes in the long range 
part of the 6 Li wave function have effects of order 1 fm 2 on the quadrupole moment. Thus 
the values of the asymptotic properties, D% d =-0.29 fm 2 and 7/ a d=-0.07±0.02, obtained with 
this wave function may not be very accurate. The total normalizations are N^ d = 0.82 and 
= 0.021. The resulting spectroscopic factor, 0.84, is in good agreement with the value 
of 0.85 obtained in radiative capture experiments p2]| . 

The two-cluster densities p^ d Id,Mj (r), multiplied by r 2 , are shown in Fig. They have 
two peaks; the smaller inner peak at r ~ 0.9 fm is almost spherically symmetric, while the 
larger peak at r ~ 4 fm is anisotropic. In particular p^ d '°(rz) is much larger than p^ d (rx) 
for r > 2 fm. In the former configuration the r is along the axis of the torus, while in the 
latter it is transverse. This anisotropy is also a consequence of the toroidal shape of the 
deuteron in the M d =0 state. 



VI. OTHER T, S CHANNELS 

In this section we discuss the properties in nuclei of pairs of nucleons with T, S — 0, 0, 
1, 0, and 1, 1. Like the T, S — 0, 1 channel discussed in the previous sections, the T, S — 1, 1 
channel also has a tensor potential, but it has the opposite sign of that for T,S = 0,1. 
Therefore the role of Ms is reversed compared to that in T = states; Ms = pairs have 
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maximum density along the z-axis, while Ms = ±1 pairs have maximum density in the 
xy plane as can be seen in Fig. |23], which shows pYf (r,9)/ Rf x for 4 He, 6 Li, and 16 0. The 



curves for Ms = ±1, 9 = n/2 are between the two sets of curves shown in the figure; to 
reduce clutter they are not shown. The curves for 4 He and 6 Li have been renormalized by 
the factors Rf x to have the same peak height as for 16 0; these factors are shown in Table II. 
We see that the shapes of the T, S = 1,1 density profiles are quite different in the different 
nuclei and that the anisotropy decreases as the number of nucleons increases. 

The analog of the deuteron in the T, S — 1, channel is the 1 So virtual bound state (VBS). 
For the Argonne Vis potential, this is a pole on the second energy sheet at E = —0.098 MeV 
or k = — 0.049i fm" 1 . Although the wave function is not normalizable, it has a local peak 
which we scale to compare to the unpolarized deuteron in Fig. ^4]. We see that it peaks at a 
slightly larger radius and is broader. The figure also shows the p? o^V-^i o of 4 He, 6 Li, and 
16 0; the curve for 3 He is between those of 4 He and 6 Li, while the curve for 7 Li is very close 
to that of 6 Li. Again the curves have been normalized to the peak height of the 16 density. 
We see that the short-range shapes of the p\ in nuclei are well reproduced by the VBS 
density. Finally, Fig. UB shows the densities for the T, S — 0, channel, again normalized to 



16 0. As is the case for T, S = 1, 1, there is no common shape. 

Table II also shows the number of pairs, N^ s , in these T, S channels and the corre- 
sponding IP values. As is the case for T, S — 0, 1 (Table I), the number of pairs increases 
more rapidly with A than does s , because of the increasing proportion of pairs with large 
separation. 

Using the projection operators (1 — Tj • Tj)/A and (3 + r j • Tj)/A for T = and 1 pairs 
we find that the total number of T = and 1 pairs in a nucleus depends only on its mass 
number A and isospin T A . 

Ko + N o,i = \[A 2 + 2A- 4T a (T a + 1)] , (6.1) 
Ni + iVA = I \3A 2 -6A + 4T A (T A + 1)1 . (6.2) 

o L J 

The above relations are obeyed by N^ s obtained from either the IP or correlated wave 
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functions, since, in the present study, both are eigenstates of T4. 
If the total spin, 

S A = '5l\ a 'i> ( 6 - 3 ) 

i 

were to be a good quantum number we would have similar relations, 

<o + K = I [ A 2 + 2A - 4S A (S A + 1)] , (6.4) 

Ki + N t = I [ 3A " ~ 6A + AS a(Sa + 1)] , (6.5) 

for the total number of pairs with spin and I. They are obeyed by the N^ s calculated 
for the IP states which have S A = 1, J, 0, 1, J and for 2 H, 3 He, 4 He, 6 Li, 7 Li and 16 
respectively. However, S A is not a good quantum number; tensor correlations admix states 
with larger S A in the ground state. These reduce the Nf and increase the Nf^ by the same 
amount due to Eq. (6.2). In 3 He ( 4 He) the is given by 1.5P D (3P D ), where P D is the 
fraction of L — 2, S A — | (Sa = 2) state in the nuclear ground state. 

The interaction in the T,S — 1,0 state is much more attractive than that in the T,S = 
1, 1 state. Hence the depletion of T, S — 1,0 pairs by tensor correlations reduces the binding 
energy of nuclei significantly. For example, in 4 He the T, S = 1,0 interaction gives —14.2 
MeV per pair, while the T,S = 1,1 interaction gives only -0.8 MeV per pair. Thus the 
conversion of 0.47 T=l pairs from S=0 to S=l state raises the energy of 4 He by ~ 6.3 
MeV. It should be stressed that this is a "many-body" effect absent in the two-body cluster 
approximation of either Brueckner or variational methods. The tensor interaction between 
nucleons % and j can flip their spins and convert pairs ik and/or jl from 5=0 to S—l. 



VII. CONCLUSIONS 

The main conclusions of this study of nuclear structure, as predicted by realistic models 
of nuclear forces, are: 

(i) The static part of the two-nucleon potential in T, S, Ms = 0, 1, state has a large 
angular dependence due to the tensor interaction dominated by one-pion exchange. At 
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r ~ 1 fm the difference between this potential at 9 = 7r/2 and is ~ 300 MeV in most 
models (Fig. [I]). It confines T, S, Ms = 0, 1, pairs to the small 9 region producing toroidal 
distributions. The central hole in these tori is due to the repulsive core in NN interaction. 
The maximum density in the tori is large, due to which the peak one-body density in 
deuterium exceeds 0.3 fm -3 in most models. 

(ii) The more familiar dumbbell (or cigar) shaped density distribution of the deuteron 
in Ms=±l states can be considered as that produced by a rotating torus. 

(iii) The diameter of the maximum density torus, and the thickness of the half-maximum 
density torus are predicted to be ~ 1.0 and 0.9 fm, respectively; these values are supported 
by the observed elastic electron-deuteron scattering. 

(iv) The pair distribution functions in T, S = 0, 1 states indicate that the tensor corre- 
lations have near maximal strength in all nuclei considered here for r < 2 fm. 

(v) The pair distribution functions in T, S = 0, 1 and 1,0 states in different nuclei, can 
be scaled to lie on universal surfaces for r < 2 fm. These universal surfaces are predicted by 
the density distributions of the deuteron and the 1 S virtual bound state. The scaling factor 
RAd for the T, S = 0, 1 densities provides a rigorous definition of the Levinger-Bethe quasi- 
deuteron number of the nucleus. The calculated values of Ra<i are significantly different from 
estimates based on independent-particle models, and in qualitative agreement with photon 
and pion absorption data. 

(vi) The many-body distribution functions are also predicted to be anisotropic. In par- 
ticular the anisotropies of the dp, dd, and ad distributions in 3 He, 4 He and 6 Li are strongly 
influenced by the toroidal structure of the deuteron. 

(vii) Tensor correlations convert T=l pairs of nucleons from S=0 to S=l state. This 
many-body effect reduces the binding energies of nuclei. It does not appear as if many-body 
effects reduce the magnitude of tensor correlations for the range of nuclei studied here: 2 H 
to 16 0. 

Due to the small size of this toroidal structure it may be worthwhile to attempt to un- 
derstand it from the more basic quark degrees of freedom. Within the constituent quark 
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model f54|j55| it requires a solution of the six-quark problem with a suitably chosen Hamil- 
tonian. Many attempts have been made (see Refs. |56|-|58| for example) to calculate the 
nucleon-nucleon interaction from approximate solutions of the six-quark Hamiltonian. A di- 
rect coupling of the pions to the quarks is used to obtain the tensor part of the interaction. 
The toroidal structure is presumably very sensitive to this coupling and to the tensor part 
of the quark-quark interaction in the framework of the constituent quark model. 

As is well known, toroidal structure for the ground state of the deuteron was predicted 
many years ago Hl6j|l7| using classical Skyrme field theory |15| related to QCD in the N c — > oo 
limit. In the classical limit one obtains a toroidal shape of ~ 1 fm in size and a binding 
energy of ~ 150 MeV. From Fig. [1] it is obvious that in the classical limit realistic models of 
nuclear forces would also give a deuteron binding energy in the 100 to 200 MeV range. There 



have been attempts to include quantum corrections to this theory. A recent calculation [18 



obtains an energy of -6. 18 MeV for the deuteron in this model. Ground states of the classical 



Skyrme field with baryon numbers 3 to 6 have also been studied |19|]. The baryon equidensity 
surfaces of these classical solutions are highly anisotropic. However, the nucleon equidensity 
surfaces of the J n = | + and + 3 He and 4 He must be spherically symmetric, thus a direct 
comparison is not possible. Nevertheless the anisotropic dp and dd distributions in 3 He and 
4 He may have some relation to the baryon density distributions in the Skyrme model. 
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FIGURES 

FIG. 1. The upper four lines show expectation values of for Ms=0, 9=0, and the lower 
four lines are for Mg=0, 9=ir/2 or equivalently Ms=±l, 9=0. The expectation values for Ms=±l, 
9=ir/2 (not shown) are half way in between. 

FIG. 2. The S- and D-wave deuteron wave functions for various potential models. 

FIG. 3. The top, middle and bottom four curves respectively show the deuteron density for the 
indicated values of Md and 9, obtained from various potential models. 

FIG. 4. The deuteron density p d (x',z') obtained from the Argonne vis model. The peaks are 
located at z'=0 and x'=±d/2. 

FIG. 5. The deuteron density p^ 1 (x', z') obtained from the Argonne v\% model. The peaks are 
located at x'=0 and z'=±d/2. 

FIG. 6. The surfaces having p| 1 (r')=0.24 fm~ 3 (A) and p^(r')=0.24 fm~ 3 (B). The surfaces 
are symmetric about z' axis and have r' <0.74 fm, i.e., the length of the dumbbell along z' axis 
as well as the diameter of the outer surface of the torus is 1.48 fm. Sections C and D are for 
p^ d (r')=0.08 fm -3 ; the maximum value of r' is 1.2 fm. 

FIG. 7. The square of the calculated deuteron charge form factors. 

FIG. 8. The values of deuteron T2o{q) obtained from Eq. (3.11) are shown by full line, whereas 
the dashed line gives T 2 o(q) including magnetic contributions for a 15° electron scattering angle. 

FIG. 9. The square of the deuteron magnetic form factor calculated with (full line) and without 
(dashed line) convection current term. 

FIG. 10. The transition density pt r (x'z') for elastic magnetic scattering by deuterons. The 
peaks are located at z'=0 and x'=±0.5 fm. 
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FIG. 11. The deuteron momentum distribution for selected values of Md and 6 k- 

FIG. 12. The calculated values of d(e, e'p)n cross section for the kinematics described in the 
text. Hollow and filled symbols indicate results of complete calculations without and with me- 
son-exchange currents. 

FIG. 13. Pof{r,0)/R A d for various nuclei. 

FIG. 14. Ro(rd p ) and R2(rd P ) for 3 He. The points show results of Monte Carlo calculations in 
configuration space, and the curves are smooth fits. The asymptotic Rl given by Eq. (5.5) are 
shown by dashed lines. 

FIG. 15. Density distribution of dp clusters in 3 He with Mj = ±, M d =0 and M p = \. The 
peaks are located at Xd p =0 and z dp ~ ±1 fm. 

FIG. 16. Density distribution of dp clusters in 3 He with Mj = §, M d =l and M p = -\. The 
peaks are located at Zd P =0 and x dp ~ ±1 fm 

FIG. 17. Momentum distribution of dp clusters in 3 He in Mj = \ state for momenta parallel 
and transverse to the z-axis 

FIG. 18. Ro(rdd) and i?2(^dd) for 4 He. See Fig. 15 for notation. 

FIG. 19. Density distribution of dd clusters in 4 He in parallel (0=0) and transverse (0=ir/2) 
directions. 

FIG. 20. Momentum distribution of dd clusters in 4 He in parallel (6=0) and transverse (6=tt/2) 
directions. 

FIG. 21. Ro(r a d) and R2(r ad ) for 6 Li. See Fig. 15 for notation. 
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FIG. 22. Density distribution of ad clusters in 6 Li in parallel (6=0) and transverse (9=ir/2) 
directions. 

FIG. 23. p^ I f(r,9)/Ri 1 for various nuclei. The upper three curves are for Ms=0, 9=0 while 
the lower ones are for Ms=0, 9=n/2 and equivalently Mg=±l, 9=0. 

FIG. 24. Pifi{r) I Ri q for various nuclei. The p(r) of an unpolarized deuteron, normalized to 
have the same maximum value, is shown for comparison by the dotted line. 

FIG. 25. Po,o( r ) / Rq,o f° r various nuclei. 
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TABLES 



TABLE I. The calculated values of Ra<i and other ratios. 



Nucleus 


R-Ad 


{v-k)a 
(Un)d 


a ab,A 
a ab,d 


a ab,A 
a lb,d 


N A 
iv 0,l 

IP 


*v 


3 He 


2.0 


2.1 


2.4(1) 


~2 


1.5 


1.49 


4 He 


4.7 


5.1 


4.3(6) 


~4 


3 


2.99 


6 Li 


6.3 


6.3 






5.5 


5.46 


7 Li 


7.2 


7.8 




6.5(5) 


6.75 


6.73 


16Q 


18.8 


22 


17(3) 


16(3) 


30 


30.1 



TABLE II. The calculated values of R$ s and N# s in various nuclei. 



Nucleus R? fi N£ Q R$ N a R^ 







IP 






IP 






IP 




3 He 


0.087 


1.5 


1.35 


0.0016 





0.01 


0.012 





0.14 


4 He 


0.22 


3 


2.5 


0.0085 





0.01 


0.060 





0.47 


6 Li 


0.24 


4.5 


4.0 


0.061 


0.5 


0.52 


0.104 


4.5 


4.96 


7 Li 


0.37 


6.75 


6.1 


0.118 


0.75 


0.77 


0.18 


6.75 


7.41 


16Q 


1 


30 


28.5 


1 


6 


6.05 


1 


54 


55.5 
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